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SUMMARY 

Droplet trajectories about a sphere in ideal fluid flow were cal- 
culated. From the calculated droplet trajectories the droplet impinge- 
ment characteristics of the sphere were determined. 

Impingement data and equations for determining the collection ef- 
ficiency, the area, and the distribution of impingement are presented 
in terms of dimensionless parameters. The range of flight and at- 
mospheric conditions covered in the calculations was extended con- 
siderably beyond the range covered by previously reported calculations 
for the sphere. 


INTRODUCTION 

One of the important factors in the design and evaluation of ther- 
mal icing-protection systems is the determination of the cloud-droplet 
impingement characteristics of the various components requiring protec- 
tion. The placement of instruments that are sensitive to impinging at- 
mospheric water droplets and ice accretions in the nose region of air- 
craft and missiles makes it increasingly important to provide protection 
for this area. An aircraft radar system, for example, is affected by a 
layer of ice or water distributed over the radome surface. 

A portion of a prolate ellipsoid of revolution is frequently a 
good approximation to the nose region of an aircraft or missile. For 
this reason, theoretical droplet impingement calculations on prolate 
ellipsoids of revolution have been made by several investigators. Im- 
pingement calculations for ellipsoids varying in fineness ratio from 1 
(a sphere) to 10 are reported in references 1 to 6. 

In addition to their use in determining the impingement of cloud 
droplets on aircraft components, the impingement characteristics of the 
ellipsoid family are of interest in other fields such as chemical en- 
gineering and physical meteorology. For example, reference 7 shows 
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that the growth of a rain drop due to accretion as it falls through 
small cloud droplets can he calculated by using the collection effi- 
ciency of a sphere. 

In contrast to the ellipsoids of other fineness ratios, the avail- 
able differential -analyzer calculations of droplet impingement on the 
sphere (ref. l) cover a very limited range of flight and atmospheric 
conditions. Therefore, it has been necessary to use theoretical and 
empirical extrapolation equations based on very limited calculated data 
in order to determine the impingement characteristics of the sphere at 
many of the physical conditions of interest. 

The calculated droplet impingement data for the sphere presented 
herein are a continuation of the study of droplet impingement on pro- 
late ellipsoids of revolution reported in references 3 and 4. The pur- 
pose of this report is to extend the range of published droplet impinge- 
ment, data for the sphere in ideal fluid flow. Droplet trajectories were 
computed at the NACA Lewis laboratory with the aid of an electromechani- 
cal differential analyzer. From the computed trajectories, the rate, 
distribution, and surface extent of impinging water were obtained and 
are summarized in terms of dimensionless parameters which cover a wide 
range of flight and atmospheric conditions. 


SYMBOLS 


°D 

d 

K 

R 

Re 

Re 0 

r ,z 

r 0 


The following symbols are used in this report: 
drag coefficient for droplets, dimensionless 
droplet diameter, microns 
collection efficiency, dimensionless 
inertia parameter, 1.704X10”^ Ud^/pR, dimensionless 
radius of sphere, ft 

local Reynolds number with respect to droplet, dimensionless 

free- stream Reynolds number with respect to droplet, 
4.813X10”^ dp a tl/p, dimensionless 

cylindrical coordinates, ratio to radius, dimensionless 

starting ordinate at z = - • of droplet trajectory, ratio to 
radius, dimensionless 
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r O,tan 


Sm 

U 

u 

v 

W 


W, 


'm 


W Q 


w 


P 

Po 


starting ordinate at z = - * of droplet trajectory tangent to 
surface of sphere, ratio to radius, dimensionless 

distance along surface of sphere from forward stagnation point 
to point of droplet impingement, ratio to radius, 
dimensionless 

limit of impingement zone, ratio to radius, dimensionless 
free- stream velocity, or flight speed, mph 
local air velocity, ratio to free-stream velocity 
local droplet velocity, ratio to free-stream velocity 
rate of impingement of water, Ib/hr 

total rate of impingement of water on surface of sphere, lb/hr 
local rate of impingement of water, lb/(hr) (sq ft) 
liquid- water content in cloud, g/cu m 

r 0 to 0 

local impingement efficiency, . dimensionless 


r dS * 

local impingement efficiency at stagnation point, dimensionless 
viscosity of air, slugs/(ft)(sec) 
p density of air, slugs/cu ft 

£L 

x time scale, dimensionless 

Subscripts: 

r radial component 

z axial component 


CALCULATION OF DROPLET TRAJECTORIES 

As discussed in reference 3, the dimensionless equations of motion 
of droplet trajectories about a body in axisymmetric flow are of the 
same form as those derived in reference 8 for two- dime nsional flow. 
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Using the same coordinate system (fig. l) as that used for the prolate 
ellipsoid in reference 3 , the dimensionless equations of motion for 
the z,r plane can he written 


( 1 ) 


K = 1.704X10" 12 (2) 

The Reynolds number with respect to the droplet Re is obtained in 
terms of the free- stream Reynolds number 

a dp„u 

Re n = 4. 813X10" b .. a _ 

from the relation 

(l^) 2 = K " v z )2 + (u r “ v r) 2 (4) 

The equations of motion (l) of the droplet trajectories in the 
z,r plane are solved and summarized in terms of the dimensionless param- 
eters ReQ and K. From the droplet trajectories in the z,r plane, 
the droplet impingement characteristics of the sphere can be obtained. 

The air-flow field around the body, required for the solution of 
the equations, was obtained from the classical equation for nonviscous, 
incompressible fluid flow about a sphere. The z- and r-components of 
the velocity can be expressed in the form 


1 (2z^ - r ) 

% ' 2 (z 2 + r 2 )5/ 2 

(5) 

3 zr 

^ “ 2 (z 2 + r 2 ) 5 / 2 

(6) 


(3) 


dv 
z 

dT 


C D Rc 1 


i K - v z ) ] 


24 K 


dv r 

dT 


C D Re 1 / , 

— K ' v r> 


where 


The values of the dimensionless air-velocity components u z and u r 
were calculated with the aid of digital computers for several hundred 
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points in the flow field. From the computed velocity values, input 
charts were constructed on glass cloth for the input tables of the 
differential analyzer used to calculate the trajectories. The values 
of the air velocity components u z and u_ are given in figure 2. 

Figure 2(a) gives u z as a function of z for constant values of r, 
while figure 2(b) gives u r as a function of r for constant values 
of z. 


The trajectories were calculated on the Lewis laboratory electro- 
mechanical differential analyzer; the characteristic features of the 
electromechan i cal differential analyzer are described in the appendix. 
The procedure for calculating trajectories of cloud droplets with re- 
spect to the sphere is basically the same as that described in refer- 
ences 3 and 8 for similar calculations on the mechanical differential 
analyzer or analog. 

The equations of motion of the droplet trajectories were solved 
for various values of the parameter ReQ between 0 and 8192. For each 

value of ReQ, a series of trajectories was computed for each of several 

values of l/K between 0.0033 and 45. In order that these dimension- 
less parameters have more physical significance in the following dis- 
cussions, some typical combinations of K and ReQ are presented in 

table I in terms of the radius and the velocity of the sphere, the 
droplet size, and the flight pressure altitude and temperature. 

Because the sphere is an aerodynamic ally blunt body, it disturbs 
the air-flow field for a considerable distance upstream. The solution 
of the droplet trajectories was performed, therefore, in three stages by 
separately solving the part of each trajectory within each of the three 
following regions: (l) z = -15 to -5, (2) z = -5 to -2, and (3) z = -2 
to 0. When the solution of the trajectories was divided into three 
parts, the optimum scale factors for each of the three regions could be 
used; thus, the accuracy was considerably improved over that which 
would be obtained with a single solution from z = -15 to 0. Tests 
showed that starting at z = -15 was equivalent to starting the drop- 
let at z = - 00 as the change in droplet ordinate and velocity between 
z = -® and z = -15 would have been negligibly small. The starting 
conditions at z = -15 were taken, therefore, to be those of the air 
stream at the chosen ordinate. The starting conditions at z = -5 and 
z = -2 were obtained from plots of the droplet velocities and ordi- 
nates (for the given ReQ and K) at the end points of the previous 
interval . 
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RESULTS AND DISCUSSION 

A series of droplet trajectories about the sphere was computed for 
the various combinations of the dimensionless parameters K and Reo 
The data obtained from these trajectories are summarized in figure 3, 
where the starting ordinate r Q of each trajectory is given as a func- 
tion of the point of impingement on the surface S (the distance meas- 
ured along the surface from the forward stagnation point (-1.0, 0, fig. 
l) to the point of impingement) . Because S is made dimensionless by 
dividing by the radius of the sphere, it is also equal to the angle 
(in radians) between the z-axis and the radius vector to the point of 
impingement. The dashed lines in figure 3 are the loci of the termini 
of the constant K curves. These loci were found to be the same, 
within the order of accuracy of the computations, for all values of 
ReQ, as can be seen by comparing figures 3(a) to (g) . Because droplets 

with infinitely large K would travel in straight lines in the free- 
stream direction, the l/K = 0 curve is obtained from the relation be- 
tween r and S for a circle. From the data presented in this figure, 
the rate, area, and distribution of water-droplet impingement on the 
surface of the sphere can be determined for given values of Reo 
and K. 


Total Rate of Impingement of Water 

In flight through clouds composed of droplets of uniform size, the 
total amount of water in droplet form impinging on the sphere is deter- 
mined by the amount of water contained in the volume within the surface 
formed by the tangent trajectories (fig. l) . Therefore, the total rate 
of impingement of water (lb/hr) can be determined from the relation 

W m = 0.33nwR 2 Ur^ taB (7) 


where 0.33 is a conversion factor, the flight speed U is in miles per 
hour, liquid-water content w is in grams per cubic meter, and R is 
in feet. When constants are combined, 

W m = 1.04wR 2 Ur^ tan (8) 


In this equation, Tq tan is equal to the collection efficiency 
(E^ = Tq i _ an ) defined as the ratio of the actual amount of water inter- 
cepted by the sphere to the total amount of water in droplet form con- 
tained in the volume swept out by the sphere. 
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The value of rp ; t.nn for a given combination of Reg and K can 
be obtained from figure 3 by determining the value of rg which cor- 
responds to the maximum S for the constant K curve of interest. 

The values of rg t flTI fall on the dashed termini curves of figure 3. 

In order to facilitate interpolation and extrapolation, the data are 

2 

replotted in the form of rg tnn as a function of K for various con- 
stant values of Reg in figure 4. Examination of figure 4 shows that 
r O,tan increases with increasing K but decreases with increasing Reg. 
At values of K less than 0.10, r^ tan approaches zero. It was not 

possible to determine from the differential-analyzer calculations 
whether or not a critical value of K exists (as suggested in ref. l) 
below which there is no droplet impingement on the sphere. However, 
for K less than the critical value given in reference 1 (K = 0.083) 
the collection efficiency is certainly very small and is approaching 
zero. 


The accuracy in the determination of r g^ tan ^ rom the calculated 

trajectories is about of the same order for the sphere as for the ellip- 
soid of fineness ratio 5 (ref. 3) . For values of r n ^ +fm greater than 

0.10, rg + nri was deteimined with an accuracy of the order of ,±0.005. 

For reported values of rp^n < 0.10, the accuracy of determining the 
tangent is more indefinite but appears to be within ±0.01. 


Extent of Droplet Impingement Zone 

The extent of the droplet impingement zone on the surface of the 
sphere is obtained from the tangent trajectories (fig. l) . The point 
of tangency determines the rearward limit of the impingement zone. The 
limit of impingement for a particular Reg and K condition can 

be determined from the maximum S value of the constant K curves of 
interest in figure 3. Again, to facilitate interpolation, the data are 
replotted in the form of Sm as a function of K for constant Reo 
values in figure 5. The data of this figure indicate that Sjjj in- 
creases with increasing K but decreases with increasing Reg. A com- 
parison of figures 4 and 5 shows that at values of K less than 0.1, 

^ approaches zero much more slowly than the collection efficiency 

O 

r 0 tan* collection efficiency is practically zero in most of the 

region between K = 0.01 and 0.10, while Sm is still quite large. 

This difference can be explained by the fact that for small values of 
K only droplets starting close to the stagnation line (-z-axis) will 



8 


MCA TN 3587 


impinge on the sphere; however, the trajectories of these droplets are 
such that in the neighborhood of the sphere their path is almost paral- 
lel to the surface. Thus, they may travel a sizeable distance along 
the surface of the sphere before impinging even though the amount of 
water collected on the surface is very small. 

Because of the difficulty of determining the exact point of tan- 
gency on the surface of the sphere of each tangent trajectory, the ac- 
curacy of determining Sm for each K and Re© condition is of the 
order of ±0.020. Because the data points were closely distributed 

about a common termini curve (fig. 3) for all values of Re Q and K, 
the error in the values of taken from the best-fit termini curve 

and replotted in figure 5 is of the order of ±0.01. The accuracy in 
determining the values of S for the intermediate points of impingement 
given in figure 3 was much higher because the points at which the in- 
termediate trajectories terminated on the surface of the sphere were 
much better defined. 


Distribution of Impinging Water Along Surface of Sphere 

The amount of water impinging on the surface of the sphere within 
any ring of width S 2 - can be determined if the starting or dina tes 

r Q are known for the droplets that impinge at Sj_ and S 2 . These 
data can be obtained from figure 3. 

The amount of water (lb/hr) impinging within the ring of width 
S 2 - Sj_ is given by the relation 

W = 1.04(r^ 2 - rg^wR^ (9) 


Local Rate of Impingement of Water 

The local rate of impingement of water in droplet form 
(lb/(hr) (sq. ft)) on the surface of the sphere can be determined from 
the expression 


r 0 ^0 


Wp = 0.33Uw — = 0.33Uwp 


( 10 ) 


where p is the local impingement efficiency. The values of p as a 
function of S for combinations of ReQ and K are presented in fig- 
ure 6. These curves were obtained by multiplying the slope of the 
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curves in figure 3 by the corresponding ratio r Q /r at each point. 
Because the slopes of the ro against S curves (fig. 3) in the re- 
gion between S = 0 and S = 0.05 are difficult to determine, the 
exact values of 3 between S = 0 and S = 0.05 are not known. The 
values of 8 presented in this region between S = 0 and S = 0.05 
are estimated to be accurate within ,±0.05. A plot of the local impinge- 
ment efficiency 3q at the stagnation point (s = 0) as a function of 

K is shown in figure 7. The local impingement efficiency at S = 0 
is sometimes useful as a parameter in empirical equations for calcu- 
lating local impingement efficiency (ref. 2). 


IMPINGEMENT IN CLOUDS OP N0NUNIF0RM DROPLET ST7.T? 

The data presented in figures 3 to 7 apply directly only to flights 
in clouds composed of droplets that are all uniform in size. The drop- 
lets in a cloud, however, may have a range of sizes. Theoretical cal- 
culations and experience in the Lewis icing research tunnel on bodies 
of revolution have shown that the amount of ice collected when a dis- 
tribution of droplet sizes is present in the tunnel is considerably dif- 
ferent from that which would be obtained if only droplets of the volume- 
median size were present. Therefore, if the size distribution of the 
cloud-droplets is known or can be estimated, then the data must be ac- 
cordingly modified (or weighted) before the rate, extent, and distribu- 
tion of droplet impingement on the sphere are calculated. Methods of 
weighting the impingement data for a given droplet-size distribution 
are described in detail in reference 9. 


CONCLUDING REMARKS 

Because the droplet trajectories were calculated for incompressible, 
ideal fluid flow, a question may arise as to their applicability at the 
higher subsonic velocities. Reference 10 shows that the effect of com- 
pressibility on the impingement characteristics of a circular cy lin der 
is small. Therefore, it is expected that the assumption of incompress- 
ibility does not invalidate the application of the calculated values at 
high subsonic velocities. A more serious deviation from ideal fluid 
flow is caused by the separation of the boundary layer from the surface 
of the sphere in the region between 70° and 90° from the stagnation 
point because of viscous effects. It is not known with certainty what 
effect the flow separation would have on droplet trajectories, but it 
is expected that except for large values of K (where the tangent tra- 
jectory touches within the separation region) the droplet impingement 
pattern would differ little from that calculated for the ideal fluid. 

At large values of K the effect is probably more pronounced, and 
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experimental tests would be required to establish the degree of agree- 
ment with ideal fluid calculations. In addition to any possible ef- 
fects on the impingement pattern on the windward side of the sphere, 
the low-pressure eddy region on the leeward side of the sphere caused 
by flow separation may entrain some of the smaller droplets and cause 
some impingement on the lee side of the sphere. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 28, 1955 
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APPENDIX - DESCRIPTION OF LEWIS ELECTROMECHANICAL 
DIFFERENTIAL ANALYZER 

The electromechanical differential analyzer is basically similar 
to the mechanical differential analyzer, or analog, described in ap- 
pendix A of reference 8. The electromechanical differential analyzer 
has 10 integrators, 10 variable coefficients, 17 fixed coefficients, 
and 8 adders, all of which are synchro- coupled through an electric 
patch board. The electric interconnections between the various units 
used in problem programing are made at the patch board. 

The variable coefficients are ball and disk integrators used as a 
multiplier by keeping the ball at one position throughout a solution. 

The ball is moved by hand cranking to the proper position on the disk 
before each solution is started. A counter, geared to the crank shaft, 
indicates the value of the coefficient. Five of these units can multi- 
ply by numbers from 0 to 2; the other five units can use factors up to 
6. By switching, the factors can be made either plus or minus. Each 
of the variable coefficients has an electronic-servo amplifier to in- 
crease the available torque for rotating its disk. 

The ten integrators each have two electronic- servo amplifiers, one 
for rotating the disk and one for moving the ball carriage across the 
disk. The lead screw of each integrator is geared to a counter that 
indicates the position of the ball on the disk. 

Each fixed coefficient consists of a synchronous repeater and a 
synchronous transmitter connected by the proper gears. 

The adders are synchro-differential motors, which assume a mechani- 
cal position proportional to the algebraic sum of two input signals. 

The motors are mechanically coupled to synchronous transmitters. 

There are three input 'tables for supplying graphical data (such as 
air-velocity components) and one large output table that can simul- 
taneously plot two functions against a single independent variable. 
Problem solutions are obtained from the output-table graphical plots 
or are read from mechanical counters coupled to the differential 
analyzer. 
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Figure 1. - Coordinate system for droplet- trajectory calculations about a sphere. 
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Figure 2. - Continued. Sphere flow field 
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Distance to point of impingement, S 
(b) Pree-stream Reynolds number, 16. 

Figure 3. - Starting ordinate as function of distance along surface to point of 
impingement . 
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Figure 3. - Continued. Starting ordinate as function of distance along surface 
to point of impingement. 
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(f) Free-stream Reynolds number, 4096. 


Figure 3. - Continued. Starting ordinate as function of distance along surface 
to point of impingement. 
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Figure 4. - Square of starting ordinate of tangent trajectory as function of inertia parameter. 
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Figure 5. - Limit of impingement zone on surface of sphere. 
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(b) Pree-stream Reynolds number, 16. 

• _ r 0 dr 0 

Figure 6. - Local impingement efficiency p = — -gg— • 
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(d) Free-stream Reynolds number, 512. 

r n dr 0 

Figure 6. - Continued. Local impingement efficiency 3 = 
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Figure 6. - Continued. Local impingement efficiency p = ^ ^0 











